TEST

# do multiple layers by not passing the first "leaflet()" call a data argument

library(tidyverse)
## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──
## ✓ ggplot2 3.3.5     ✓ purrr   0.3.4
## ✓ tibble  3.1.2     ✓ dplyr   1.0.7
## ✓ tidyr   1.1.3     ✓ stringr 1.4.0
## ✓ readr   1.4.0     ✓ forcats 0.5.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(leaflet)
parks <- sf::read_sf("../data/State_Park_Locations.shp")
streams <- sf::read_sf("../data/Streams_303_d_.shp")

leaflet() %>%
  setView(lng = -96.703090, lat = 40.81928, zoom = 10) %>% 
  addProviderTiles(providers$Esri.NatGeoWorldMap) %>%
  addMarkers(data = parks, popup = ~AreaName, label = ~AreaName) %>% 
  addPolylines(data = streams, color = "blue", 
               popup = ~paste0(Waterbody_, " - ", Impairment))

TEST